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Abstract 

The next-to-leading order and leading order analysis are performed on the 
differential cross section ratio from Drell-Yan process. It is found that the effect of 
next-to-leading order corrections can be negligible on the differential cross section 
ratios as a function of the quark momentum fraction in the beam proton and the 
target nuclei for the current Fermilab and future lower beam proton energy. The 
nuclear Drell-Yan reaction is an ideal tool to study the energy loss of the fast 
quark moving through cold nuclei. In the leading order analysis, the theoretical 
results with quark energy loss are in good agreement with the Fermilab E866 
experimental data on the Drell-Yan differential cross section ratios as a function 
of the momentum fraction of the target parton. It is shown that the quark energy 
loss effect has significant impact on the Drell-Yan differential cross section ratios. 
The nuclear Drell-Yan experiment at current Fermilab and future lower energy 
proton beam can not provide us with more information on the nuclear sea quark 
distribution. 
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1 Introduction 

In 1983, the European Muon Collaboration(EMC) reported that the ratios of 
nuclear structure functions for the iron and deuteron nuclei is not equal to unity, but a 
function of the Bjorken scaling variable x in charged lepton-nucleus deep inelastic scat- 
tering. It shows that the parton distributions are modified in the nuclear environment 
because the structure function describes the quark momentum distributions in bound 
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nucleons. From then on, the quark and gluon distributions in hadrons and nuclei have 
been one of the most active frontiers in nuclear physics and particle physics. The nu- 
clear parton distribution directly affects the interpretation of the data collected from 
the nuclear reactions at high energies'^l. The precise nuclear parton distribution func- 
tion should be very important in finding the new physics phenomena and determining 
the electro-weak parameters, neutrino masses and mixing angles in neutrino physics. 

In early 1970, S.Drell and T.M.Yan first studied the production of large-mass lep- 
ton pairs from hadron-hadron inelastic collisions, which is so-called Drell-Yan process'^' . 
According to the parton model, the process is induced by the annihilation of a quark- 
antiquark pair into a virtual photon which subsequently decays into a lepton pair. The 
nuclear Drell-Yan process of proton-nucleus collisions therefore is closely related to the 
quark distribution functions in nuclei. It is further natural to expect that the Drell-Yan 
reaction is a complementary tool to probe the structure of hadron and nuclei. 

In 1990, Fermilab Experiment772(E772) Wmade the measurement of the nuclear 
dependence of the Drell-Yan process by using SOOGeV protons bombarding on D, C, 
Ca, Fe, and W nuclei. Muon pairs were recorded in the range AGeV < M < 9GeV 
and M > llGeV. The covered kinematical ranges were 0.1 < x < 0.3, where x is the 
momentum fraction of the target parton. Their aim is to investigate the modification 
of the quark structure in nuclei. The theoretical models, which can well describe the 
EMC effect in the charged lepton-nuclei deep inelastic scattering, were used to fit the 
observations on the Drell-Yan differential cross section ratios. It is found that some of 
the theoretical models overestimate the nuclear Drell-Yan ratios from Fermilab E772, 
such as pion-excess model '^1 and quark-cluster model '^1 . It indicated that whether there 
is another nuclear effect apart from the nuclear effects on the parton distributions as 
in charged lepton-nucleus deep inelastic scattering. 

In 1999, Fermilab Experiment866(E866) '^1 performed the precise observation 
of the ratios of the Drell-Yan cross section per nucleon for an SOOGeV proton beam 
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incident on Be, Fe and W target. The Drell-Yan events extend over the ranges AGeV < 
M < SAGeV, 0.01 < X2 < 0.12, 0.21 < xi < 0.95 and 0.13 < xp < 0.93, where Xi(2) 
and xf are the momentum fraction of the beam parton (the target nuclear parton) 
and Feynman scahng variable, respectively. The extended kinematic coverage of E866 
significantly increases its sensitivity to the nuclear shadowing effect and the quark 
energy loss. This is the first experiment on the energy loss of the projectile particle 
moving through cold nuclei. The E866 compared their experimental data with the 
results from leading-order Drell-Yan calculations by using the EKRS nuclear parton 
distributions'^! together with the MRST parton distribution functions'^!. It is shown 
that the energy loss effect can be negligible. 

In previous works'^^^l, we investigated the Drell-Yan differential cross section 
ratios as the function of the momentum fraction of the beam parton from E866 data 
in the framework of Glauber model by means of EKRS and HKM nuclear parton dis- 
tribution functions'^^l . We found that the theoretical results with energy loss of the 
beam proton are in good agreement with the Fermilab E866 experiment by means of 
HKM nuclear parton distributions. However, the calculated results without energy loss 
can give good fits by using EKRS nuclear parton distribution functions. Furthermore, 
we introduced two typical kinds of quark energy loss parametrization, i.e.the linear 
and quadratic quark energy loss with the average path length of the incident quark 
in the nucleus A. The nuclear dependence of the Drell-Yan production cross sections 
were calculated by combining the quark energy loss effect with the EKRS, HKM and 
HKN'^^! nuclear parton distribution. The global fit to the Drell-Yan differential 
cross section ratios indicated that the theoretical results without energy loss effects 
agree very well with the E866 experimental data by taking advantage of the EKRS nu- 
clear parton distribution functions. If employing the HKM nuclear parton distribution 
function, the results with energy loss effect are in good agreement with the Fermilab 
E866 data. In addition, the results with HKM are most near to those with HKN. It 
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is demonstrated that at current Fermilab incident proton energy, we can not distin- 
guish between the hnear and quadratic dependence of the quark energy loss. Further 
experiments are needed about nuclear Drell-Yan reactions with a lower energy incident 
proton. Using the values of quark energy loss from a fit to E866 experimental data, 
the prospects are given for the lower energy proton beam bombarding deutcron and 
tungsten. It is shown that these future experiments can give valuable insight in the 
enengy loss of fast quark propagating through cold nuclei. In the analysis above, we 
performed the leading order calculations of the Drell-Yan differential cross sections. 
As is well known, QCD corrections can alter quite significantly the cross sections at a 
hadronic collider. Thus, these may have serious bearing on the discovery potential of 
the Drell-Yan reactions, in which the leading order (LO) results may seriously underes- 
timate the cross sections. This has led to the incorporation of the next-to-leading order 
(NLO)results. As for the nuclear parton distributions, EKRS, HKM and HKN nuclear 
parton distributions are obtained by Eskola et al.'^l, Hirai et al.'^^'^^l with the global 
analysis of the relative experimental data, respectively. It is noticeable that HKM use 
only the charged lepton-nucleus deep inelastic scattering experimental data, HKN and 
EKRS employ E772 and E866 nuclear Drell-Yan reaction data in order to pin down 
the nuclear antiquark distribution in the small x region. In this report, we will explore 
the effect of NLO correction on the Drell-Yan differential cross section ratios, and the 
influence of quark energy loss on the Drell-Yan differential cross section ratios as a 
function of the momentum fraction of the target parton. 

This paper is organized as follows. In sect. 2, a brief formulism for the differential 
cross section in the Drell-Yan process is presented. The effect of NLO correction is 
given on the Drell-Yan differential cross section ratios in sect. 3. The influence of quark 
energy loss is given on the Drell-Yan differential cross section ratios as a function of 
the momentum fraction of the target parton in sect. 4, and the summary is given in 
sect. 5. 
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2 Brief formulism for differential cross section in Drell-Yan reaction 

In the Drell-Yan process'^l , the perturbative QCD leading order contribution is 
quark- ant iquark annihilation into a lepton pair. The annihilation cross section can be 
obtained from the e+e~ — > h'^/jl' cross section by including the color factor | with the 

charge for the quark of flavor /. 

da Svra^ n c-/- , /.x 

7m = ^'''^'-''>' 

where V§ — (xiX2s)^/^, is the center of mass system (cm. system) energy of qq collision, 
a;i(resp.a;2)is the momentum fraction carried by the projectile (resp. target) parton, ^/s 
is the center of mass energy of the hadronic collision, and M is the invariant mass of the 
produced dimuon. The hadronic Drell-Yan differential cross section is then obtained 
from the convolution of the above partonic cross section with the quark distributions 
in the beam and target, 

where a is the fine-structure constant, the sum is carried out over the light fiavor 
f = u,d, s, and q^^'^\x, Q'^) and (ff^^\x, Q'^) are the quark and anti-quark distributions 
in the proton (nucleon in the nucleus A). 

In addition to the leading order Drell-Yan term, the contributions are needed 
from qq annihilation processes( qq ^ J* + g) and gluon Compton scattering( q-\- g ^ 
7* + ?)'^^' in the perturbative QCD next-to-leading order, which are denoted as Ann 
and Comp, respectively. The contribution from the order-a^ annihilation graphs is 

j2 Ann „i J2-Ann 

j2 ^ Ann 

+^^{ti--h)^e}qj{t,,Q')qf{h,Q% (3) 



with 
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Sa^a.iQ-) + 1 2,ri(l-,r2) 

+ ^ 27Q- - kih-xih X2{ti+xi) 

3 2 3xi, ^, 

+257^-^-^^ + ^^-^) 

^ ^ 16«2^,(Q2) ^ (r + t,t2)[r2 + (tits)'] 



f 27g2 ^(tit2)2(ti+Xi)(t2 + X2)[(tl-Xi)(t2-X2) 

2r(r + tits) 1 



(4) 



tit2(tiX2 + t2Xi)2 

The contribution from the Compton scattering graphs is 

j2 Comp „i j2±.C'omp 

j2±.Comp 

- ^2) E ej[g^(ti, g2) + g^(ti, g')]5^(t2, Q')}, (5) 

where g{t, Q^)is the gluon distribution in the beam proton and target nuclei, and 
rf'^^°"^^ ^ 2a2a^ _ x^ + (t2-xi)^ 2xi(l-X2) 

dtldt2 y ^^^2 ^2j[ X2(tl+Xi) 

J- 3xi(ti - Xl) 

^2ti tf J 

2a2a,(Q2) ^.^(r + tit2)[r2 + (r - tit2)^] 



+ E 



Y 9(52 tfti(tia;2 + t2a:^i)(t2 + a:;2)(t2 - 3:^2)+ 

t(t + tita) [tit|xi + T(ti3;2 + 2t2a:i)] 
(tit2)2(tiX2 + t2a;i)3 J- 

Therefore, to the next leading order, the differential cross section in Drell-Yan reaction 
can be written as 

dxidx2 dxidx2 dxidx2 dxidx2 
With calculating the integral of the differential cross section above in leading order 
(Eq.2) and next-to-leading order (Eq.7), the Drell-Yan production cross section is given 

by 

"^"^ ^ d^2{i)-^^4:r- (8) 



dxn2) J dx\dx2 



3 The influence of QCD correction on the ratio of DreU-Yan cross section 
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In the Drell-Yan reaction experiments, the ratios are measured of Drell-Yan cross 
sections on two different nuclear targets bombarded by proton, 

In order to explore the influence of QCD correction, we introduce the ratios, 

KfirM = <.%(^i(2))/<%.(^i(2)), (10) 

wherei?^^^/^2(a;i(2)) and -R;ai%2(^i(2)) ratios of Drell-Yan differential cross sec- 

tion with QCD correction and with only leading order contribution, respectively 

By taking advantage of the HKM cubic type of nuclear parton distribution 
functions ^-^^^ the ratios -R^j^*^{^^(a;i(2)) and -R^/^g^'^(a;i(2)) for proton incident on 
Be, Fe, W targets are calculated at the Fermilab SOOGeV proton beam energy in the 
range 4 < M < 8GeV. It is found that the differential cross section ratios in the 
next-to-leading order are almost identical to those in the leading order. The similar 
results are given for the lower energy proton bombarding deuterium and tungsten 
at the Fermilab Main Injector (FMl,120GeV proton beam)!^"^' and the Japan Hadron 
Facility(JHF,50GeV proton beam)[^^l Therefore, it can be concluded that the QCD 
correction is negligible in the nuclear Drell-Yan reactions for the current Fermilab and 
lower energy proton beam . The production of lepton pairs in proton-nucleus collisions, 
the nuclear Drell-Yan process, is one of most powerful tools to probe the propagating of 
quark through cold nuclei. The experimental study of the relatively low energy nuclear 
Drell-Yan process can give valuable insight in the quark energy loss dependence on the 
medium size. Furthermore, the influence of quark energy loss can be investigated on 
the Drell-Yan differential cross section ratios as a function of the momentum fraction 
of the target parton in the leading order. 



4 The impact of quark energy loss on the Drell-Yan diflFerential cross 
section ratios as a function of the momentum fraction of the target peirton 
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The Fermilab E866 collaboration reported their measurement of the differential 
cross section ratios Rpe/Bei^T) and Rw/Bei^-i) for an SOOGeV proton beam bombarding 
Be, Fe and W nuclei '^1 . By combining HKM cubic type of nuclear parton distribution!^^! 
with the quark energy loss, the global analysis to the E866 experimental data are 
performed in the perturbativc QCD leading order. We introduce two typical kinds of 
quark energy loss expressions. One is rewritten as 

Ax, = a^^, (11) 

where a denote the incident quark energy loss per unit length in nuclear matter, < L >a 
is the average path length of the incident quark in the nucleus A, Ep is the energy of 
the incident proton. The average path length is employed using the conventional value, 
< L >A= 3/4(1. 274^/^^fm. In addition to the hnear quark energy loss rate, another one 
is expressed as 

Ax:=/3^4^, (12) 

i.e., the quark energy loss is quadratic with the path length. With considering the 
quark energy loss in nuclei, the incident quark momentum fraction can be shifted from 
x'l = 2^1 + Axi to Xi at the point of fusion. 

If the quark energy loss do not put in, the obtained value is = 48.06 for 
the 16 total data points. The per degrees of freedom is given by /d.o.f. = 3.00. 
It is apparent that theoretical results without energy loss effects do not significantly 
agree with the E866 experimental data on the ratios Rai/A2{.^2)- After adding the 
fast quark energy loss, the per degrees of freedom is ^ jd.o.j . — 1.12 for the linear 
quark energy loss formula with a — 1.27, and the per degrees of freedom is given 
by /d.o.f. = 1.13 for the quadratic quark energy loss expression with /3 = 0.19. By 
using the parameter values in the linear and quadratic expressions obtained by fitting 
the E866 experimental data on the ratios -Rai/^2(^i)' P^r degrees of freedom are 

X^/d.o.f. = 1.69 with a = 1.99 and x^/d-O-f- = 1-65 with /3 = 0.29, respectively. The 



calculated results with the linear and quadratic energy loss expression are shown in 
Fig.l and Fig.2. which is the Drell-Yan cross section ratios for Fe to Be and W to Be as 
functions of X2, respectively. The solid curves are the ratios with only the nuclear effect 
on the parton distribution as in deep inelastic scattering, the dotted and dash curves 
respectively correspond to the linear energy loss with a = 1.27 and a = 1.99 in Fig.l, 
and the quadratic energy loss with P — 0.19 and P — 0.29 in Fig.2. From comparison 
with the experimental data, it is shown that our theoretical results with energy loss 
effect are in good agreement with the Fermilab E866. Meanwhile, It is noticeable that 
the values of the parameter a {or P) are different by means of the global analysis to 
the E866 experimental data on the ratios Raj^/A2{xi) and Raj^/A2{x2)- The results may 
be originated from the experimental precision. If the experimental data are sufficiently 
precise, the values of the parameter a (or /3) in the quark energy loss expression should 
be the same for fitting the ratios i?Ai/A2(^i) or Rai/A2{x2) from the nuclear Drell-Yan 
experiment. 

In order to clarify the impact of quark energy loss on the Drell-Yan differential 
cross section ratios as a function of the momentum fraction of the target parton, the 
ratios on Rai/A2{^2) without quark energy loss to those with linear quark energy loss 
are calculated and tabulated in Table 1. for the kinematic ranges covered by the E866 
experiment. The similar results can be obtained for the quadratic quark energy loss. 
It is shown that the quark energy loss effect has obvious influence on the Drell-Yan 
differential cross section ratios Rai/A2{^2)- As for the ratios RFe/Be{x2), the variation 
is approximately 2% to 4%. The extent to which the ratios Rw/Be{x2) vary is roughly 
4% to 8%. It can be deduced that the ratios of Drell-Yan differential cross section for 
nuclei A versus deuterium, Ra/d{x2), are affected largely because of the quark energy 
loss effect. In the global analysis of nuclear parton distribution functions, the Drell-Yan 
data are taken mainly to determine the sea quark modification in the small x region. 
It is obvious that, considering the existence of quark energy loss, the application of 
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TABLE 1: the ratios of RA^/A^i'^i) without quark energy loss to those with hnear 
quark energy loss. 



X2 




0.03 


0.05 


0.07 


0.09 


0.12 


a(1.27) 


Fe/Be 
W/Be 


1.018 
1.038 


1.017 
1.036 


1.019 
1.040 


1.022 
1.045 


1.027 
1.055 


q;(1.99) 


Fe/Be 
W/Be 


1.029 
1.061 


1.027 
1.058 


1.030 
1.064 


1.035 
1.072 


1.042 
1.088 



nuclear Drell-Yan data is remarkably subject to difficulty in the constraints of the 
nuclear antiquark distribution. 

The ratios of Drell-Yan differential cross section for tungsten versus deuterium, 
Rw/d{x2)i are also calculated at 50GeV and 120GeV proton beam by means of the 
HKM cubic type of nuclear parton distribution functions '^^1 and two kind of quark 
energy loss expressions. It is indicated that the theoretical results with quark energy 
loss effect deviate significantly those with only the nuclear effects on the structure 
function. As an example, Figure 3 shows the ratios Rw/d{x2) with and without linear 
quark energy loss at 50GeV and 120GeV proton beam, respectively. The kinematic 
ranges cover 4 < M < '&GeV in order to avoid contamination from charmonium de- 
cays. In these calculation, the energy loss per unit length a equal 1.27GeV/ fm and 
1.99GeV/ fm from a good fit to E866 data. Therefore, the Drell-Yan experiment at a 
lower energy proton beam do not provide us with more information on the nuclear sea 
quark distribution, especially for large-x region. 

In addition, we calculate the ratios Rw/d{^'2) with the linear and quadratic 
energy loss formula in the 50GeV and 120GeV proton beam bombarding deuterium 
and tungsten target. It is presented that the results with linear quark energy loss 
are almost identical to those with quadratic energy loss. The ratios Rw/d{x2) do 
not determine whether the energy loss is linear or quadratic with the path length. 
However, The ratios Rw/ nixi) can easily distinguish between L and dependence 
of quark energy loss^^^l. As an example, the Fig.4 shows the ratios Rw/d{x2) of the 
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Drell-Yan cross section per nucleon for an 120GeV proton beam incident on D and 
W target, where the sohd and dotted hnes correspond to a quadratic energy loss of 
P = 0.19GeV/ fm^ and to a Unear energy loss of o; = 1.27GeV/ fm, respectively 

5 Concluding remarks 

As a summary, the next-to-leading order and leading order analysis are per- 
formed on the differential cross section ratios from nuclear Drell-Yan process. The 
calculated results indicated that the QCD correction can be ignored in the nuclear 
Drell-Yan reactions for the current Fermilab and lower energy proton beam. Based 
on this view, the nuclear Drell-Yan process is one of most powerful tools to probe the 
propagating of quark through cold nuclei. The experimental study in the relatively 
low energy nuclear Drell-Yan process can give valuable insight in the quark energy loss 
dependence on the medium size. Furthermore, we have made a leading-order anal- 
ysis of E866 data on the Drell-Yan differential cross section ratios as a function of 
the momentum fraction of the target parton by taking into account of the energy loss 
effect of fast quarks. It is found that the theoretical results with quark energy loss 
are in good agreement with the Fermilab E866 experiment. The quark energy loss 
effect has obvious impact on the Drell-Yan differential cross section ratios RAi/A2i^2)- 
The nuclear Drell-Yan experiment at current Fermilab and future lower energy proton 
beam can not provide us with more information on the extraction of nuclear sea quark 
distribution. In fact, by means of the structure function xF^i^XjQ^) in neutrino deep 
inelastic scattering only, the nuclear modifications to the valence quark distribution can 
very precisely be determined in the medium and large x regions. With the structure 
functions ^2(0;, Q^) from the neutrino and charged- leptons deep inelastic scattering off 
nuclei, the nuclear modifications to the sea quark distribution in the medium and large 
X regions would be pinned down^^^l. The nuclear sea quark distribution in the small-a; 
regions can be relatively determined from the constraints of momentum conservation 
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and the structure function in small-x range. We desire to operate precise measurements 
of the experimental study in the relatively low energy nuclear Drell-Yan process. These 
new experimental data can shed light on the energy loss of fast quark propagating in 
a cold nuclei. 
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Figure caption 

Fig.l The nuclear Drell-Yan cross section ratios Rai/A2{^2) on Fe to Be(left) and W 
to Be (right). Sohd curves correspond to nuclear effects on structure function . Dotted 
and dash curves show the combination of HKM cubic type of nuclear parton distribu- 
tions with the linear energy loss a = 1.27 and 1.99, respectively. The experimental 
data are taken from the E866[7]. 

Fig. 2 The nuclear Drell-Yan cross section ratios RAi/A2i^2) on Fe to Be(left) and 
W to Be(right). Solid curves correspond to nuclear effects on structure function . 
Dotted and dash curves show the combination of HKM cubic type of nuclear parton 
distributions with the quadratic energy loss P — 0.19 and 0.29, respectively. The 
experimental data are taken from the E866[7]. 

Fig. 3 The nuclear Drell-Yan cross section ratios -^^1/^2(3^2) on W to D at 120GeV 
and 50GeV incident proton beam with the linear energy loss a = 1.27GeV/ fm{dash 
curves) and a — 1.99Gey//m (dotted curves). Sohd curves correspond to nuclear 
effects on the structure function. 

Fig. 4 The nuclear Drell-Yan cross section ratios Ra-^/A2{^2) on W to D at 120GeV 
incident proton beam with the quadratic energy loss/3 = 0.19Gey^//m^ (solid curve) 
and the linear energy loss a — 1.27GeV/ fm{dotted curve). 
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